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Abstract. - We consider the full driven quantum dynamics of a qubit realized as spin of electron in 
a one-dimensional double quantum dot with spin-orbit coupling. The driving perturbation is taken 
in the form of a single half-period pulse of electric field. Spin-orbit coupling leads to a nontrivial 
evolution in the spin and charge densities making the dynamics in both quantities irregular. As a 
result, the charge density distribution becomes strongly spin-dependent. The transition from the 
field-induced tunneling to the strong coupling regime is clearly seen in the charge and spin channels. 
These results can be important for the understanding of the techniques for the spin manipulation 
in nanostructures. 
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Introduction. Driven qubits (quantum two-level systems) attract attention due to the 
richness of phenomena occurring in different regimes of coupling to external field [T] and 
possible applications for quantum information devices. Spins of electrons in semiconductor 
quantum dots (QDs) are widely expected as a possible realization of qubits. [2] In the 
presence of spin-orbit (SO) coupling even simple systems such as the single-electron QDs 
show a rich dynamics in coupled charge and spin channels. The spin dynamics in relatively 
weak electric fields in a single QD is well understood in terms of the electric dipole spin 
resonance (EDSR). In this process the electric field at the frequency matching the Zeeman 
resonance for electron spin in magnetic field drives the orbital motion, and as a result, 
causes a spin-flip. This coupling is in the core of the proposal by Rashba and Efros [3l|4] 
for a new technique to manipulate the spin states by an external electric field. The EDSR 
spin-flip rate is orders of magnitude greater than the rate of the transitions due to the 
magnetic component of electromagnetic field. The efficiency of EDSR for GaAs QDs has 
been proven experimentally in Ref. [5] where the electric-field induced spin Rabi oscillations 
have been observed, demonstrating the abilities of a coherent spin manipulation. Another 
approach was employed in Ref. [6], where the electric field was used for driving electrons 
in a nonuniform magnetic field, thus causing a spin dynamics. The coherent spin responce 
shows that the problem with the spin dephasing in GaAs QDs, also arising due to the SO 
coupling, [?1|S] can be overcome. 

A generic system to study the driven behavior is a single parabolic one-dimensional 
QD where the electron states are the states of the harmonic oscillator. The coordinate 
and momentum can be expressed with the ladder operators satisfying known commutation 
relations. As a result, well-known selection rules for coordinate and momentum matrix 
elements can be employed and the consideration can be sufficiently simplified. However, 
completely isolated QDs cannot be used for quantum information applications since an 
intcrdot interaction is necessary to produce and manipulate many-body states and a large 
scale of the system is required to make it work. Moreover, the experiments necessarily use 
at least a double QD [36] to detect the driven spin state relative to the spin of the reference 
electron, and the role of the induced interdot tunneling on the spin Rabi oscillations has 
been noticed. [5] This "photon-induced" tunneling can put a limitation on the abilities to 
manipulate the spins. 

The information on the coupled quantum dynamics of orbital and spin degrees of freedom 
when both are nontrivial, is, however, scarce. Here we make a step forward and address full 
driven spin/charge quantum dynamics in a one-dimensional double QD |WU2| by considering 
a pumped motion of electron. These systems lavished attention due to their deceptive 
simplicity, where despite a well-established Hamiltonian, the rich variety of phenomena can 
be observed. In the classical regime, where only over-the-barrier motion can occur, the 
dynamics was studied in Ref. [13j and revealed interesting coupled irregular behavior of 
the spin and charge motion. In the more experimentally interesting quantum regime, the 
tunneling between single QDs plays the crucial role in the low-energy states, while the 
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Fig. 1: Free evolution of cr}{(t) for all spin components, as marked near the plots for the initial 
spin-up case: (a) for weak field H~ = 0.1 T (here we use a = 1.0 ■ 10~ 9 eVcm, /3 = 0), (b)-relatively 
strong field H z = 2 T. The probability to find the electron in the right dot wn(t) almost coincides 
with crfj(i), and only a slight spin deviation from the initial value due to SO coupling is observed. 

higher-energy ones are extended at a larger spatial scale. In the presence of SO coupling the 
electron states and the interdot tunneling become spin-dependent. For the driven systems 
this spin-dependence will be seen below. Nonlinear systems driven by a time-dependent 
field can show a strongly irregular chaos-like behavior. [T^HTTj] However, the existence of 
the quantum chaos, in contrast to the classical one, is disputed due to a finite set of energy 
cigenstatcs of the system. At the same time, the interplay of irregularities for spin and 
charge degrees of freedom is important for the entire system dynamics: the driven motion 
being not chaotic, is, however, strongly irregular. These irregularities, both in the orbital 
and the spin motion, limiting the abilities of the spin manipulations, will be also of our 
interest in this paper. 

Hamiltonian and time evolution. We describe a one-dimensional double QD by the 
Hamiltonian, similar to suggested in Ref. j!7) . for the particle with mass m in a quartic 
potential 



where the minima located at d and — d are separated by a barrier of height Uq. We assume 
that the barrier is opaque enough to ensure that the ground state is close to the even linear 
combination of the oscillator states located near the minima. The frequency of the single-dot 
harmonic oscillator is determined by ujq = 2\J2U^j d\fm and the corresponding oscillator 



U{x) = U (-2(x/d) 2 + {x/df) 



(1) 
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Fig. 2: Evolution of the spin components (as marked near the plots) for initial spin-up state in a 
weak magnetic field H = 0.1 T, where the spin is influenced only by the SO coupling. As in figdja), 
here we use a = 1.0- 10 -9 eVcm, /3 = 0. The external field is weak, / = 1/8, and the pulse duration 
T — 25 ps. 



length ciq = y/h/mojQ. The intcrdot semiclassical tunneling probability in the ground state: 

8^2 y/mU^d 



Wt = cxp 



(2) 



is small. The tunneling splits the ground state with the energy close to — J7o + by 
AE g of the order of HuiQWt- The high-energy states with the orbital quantum number n ^> 1 
can be treated semiclassically for the potential U{x) = Uo(x/d) 4 , yielding the eigenenergy: 

where T is the gamma function. Correspondingly, the maximum semiclassical momentum 



achieved when the electron passes the vicinity of the x — point, is p™ ax ~ \/2mE n ~ 
n 2 / 3 (hw/U ) 1/6 h/oo. 

For the system in a static magnetic field H z and driven by an external electric field £(t) 
directed along the a;-axis, the Hamiltonian is the sum of three terms H = Hq + H so + V 
with: 

2 

#o = ^ + U{x) + 9 -^ B a z H z (4) 
H so = (j3a x + a<jy)p x , V = -e£(t)x, (5) 

where e is the electron charge, and Oi are the Pauli matrices. The effect of H z is described 
by the Zeeman term only, where g is the Lande factor. The SO interaction is the sum of 
bulk-originated Dresselhaus (/3) and structure-related Rashba (a) terms. The corresponding 
velocity operator 

v = x = + H so , x] = p x /m + (3a x + uo y . (6) 

is spin-dependent. 

To describe the driven motion we use the following approach. First, we diagonalize ex- 
actly the Hamiltonian Hq+H so in the truncated basis of spinors ipn{x) |f) with corresponding 



p-4 



Pulsc-pumpcd double quantum dot 



eigenvalues E na and obtain the new basis set \ip a ) where bold n incorporates the corre- 
sponding spin index. The first four low-energy states are then defined as: |V>i) = ipi(x) |t), 

1^2) = M*) 14), 1-03) = Ms) It), and IV4) = Mm) It)- 

Since the energy E ns increases as n 4 / 3 , the truncated basis is sufficient for the purposes 
we consider below. Then, we build in this basis the matrix of the Hamiltonian V and study 
the full dynamics with the wavcfunctions in the form: 

\*)=J2Ut)e- iE » t/H W n ). (7) 

II 

The time dependence of £ n (t) with V = —e£(t)x is calculated as: 

=iZ£(t)^U(t)x«*>e-* E -- E -W\ (8) 

where x nm = (tp a \ x \ipm)- Since we do not assume a periodic driving field, we do not resort 
to the Floquet method [T2"lll8j . but use the direct numerical calculation instead. Taking into 
account that we are interested in a relatively short-term dynamics, we neglect momentum 
relaxation. This approximation is allowed at low temperatures and low excitation energies. 

The spin-dependence of the matrix element of coordinate responsible for the spin dynam- 
ics can be seen from expression i(E n — E m )x nm = H{ip n \v\ip m ) , and the spin-dependent 
terms in the velocity in ([6]). 

Observables and results. Our main goal is to calculate the probability- 

p(x,t)=tf(x,t)V(x,t) (9) 

and spin- 

Si(x,t) = &(x,t)ai$>(x,t), (10) 

density using J5]). With these distributions we find the gross quantities for the right QD: 

f 00 

u>B.(t) = / p(x,t)dx, (11) 
Jo 

poo 

oW) = / Si(x,t)dx, (12) 
Jo 

where wn(t) is the probability to find electron and <J R (t) is the expectation value of spin 
component. As the electron wavefunction at t = we take linear combinations of four 
low-energy states. The initial state is localized in the left QD, and we choose combinations: 
£1 = £3 = 1/V2 for the spin-up, and £2 = £4 = 1/V2 for the spin-down states. 

We consider a structure with d = 25\/2 nm and Uq = 10 meV, with the corresponding 
oscillator length ao = 23 nm. The four lowest spin-degenerate energy levels are E\ = 3.938 
meV, E 2 = 4.030 meV, E 3 = 9.782 meV, £ 4 =11.590 meV counted from the bottom of 
a single QD. The tunneling splitting AE g = E2 — E\ = 0.092 meV. The parameters of 
the SO coupling are assumed to be a = 1.0 • 10~ 9 eVcm, /3 — 0.3 • 10~ 9 eVcm. We 
assume the magnetic field H z = 2 T, which produces the Zecman spin splitting of the levels 
A z = E„i — E n f = 0.054 meV for the Lande-factor of electron in GaAs g = —0.45. 
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Fig. 3: Evolution of the components of the wave function (as marked near the plots) for initial 
spin-down state. The external field is characterized by / = 1/2, and the pulse duration T = 25 ps. 



In numerical calculations the basis of 20 states with the energies up to 42 meV was 
employed. We consider the period of time evolution < t < 100 ps tracking the system in 
the transient and following stationary processes. 

We begin with the analysis of the free evolution, where V — 0. The results are presented 
in figH] for two magnetic fields. One can see the free tunneling Rabi oscillations with the 
period Tr of approximately 50 ps and the corresponding spin oscillations. Figure [T] shows 
that no spin flip can be achieved in the tunneling at this system parameters. Moreover, wji(t) 
(not shown in the figure) and erfj(t) almost coincide, confirming that tunneling by itself docs 
not lead to a considerable spin dynamics. The similarity of the free motion pictures for 
H z =0.1 and H z = 2 T shows that the z— axis field, even producing the Zeeman splitting 
comparable with the tunneling splitting of the ground state, does not modify the tunneling. 
In the tunneling, spin precession mainly around the y-axis due to the Rashba SO term 
is observed. The precession angle defined for this case as <f> = atan (c^/cfj) achieves the 
maximum (f> max ~ 0.25 at t = Tr/2. This number should be compared with the precession 
angle for the classical motion, where displacement of electron by the distance I leads to 
the precession angle <f> c i(l) = 2mla/h 2 . For given geometry of the structure with I = 2d, 
one obtains <p c i(2d) — 1.2, considerably larger than max - This implies that combination of 
tunneling and SO coupling cannot be treated semiclassically being strongly different from 
what can be expected. Namely, due to the enhanced role of the tunneling the precession 
angle becomes considerably smaller. Here a general remark is appropriate. It was suggested 
[T§] and analyzed J201I1I] how to use the Larmor spin precession in magnetic field as the 
measure of the time electron spends inside a barrier. In the presence of SO coupling every 
electron has an associated momentum-dependent magnetic field acting on the spin. Whether 
the spin precession in this field can be used as a measure of tunneling time is an interesting 
question. 

As the external perturbation we consider half-period electric field pulses in the form 
£{t) = £ sm(irt/T) for < t < T. For the pulse duration we assume T = Tr/2. The 
spectral width of the pulse covers both the spin and the tunneling splitting of the ground 
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Fig. 4: Evolution of the probabilities wn(t) (a) and cr^t) (b) under the electric field pulse with 
/ = 1/8, / = 1/2, and / = 2 for the initial spin- up state. The pulse durations T = 25 ps. Dashed 
line shows the / = 1/8 regime. With the increase in the pulse field, the motion becomes less regular, 
and more high-energy component contribute into the dynamics. 

states, thus, driving the spin and orbital dynamics simultaneously. Chosen pulse duration 
satisfies the condition Tujq 3> 1 , that is the corresponding frequencies are much less than the 
energy difference between the orbital levels corresponding to a single dot, and, therefore, the 
high energy states follow the perturbation adiabatically. The field strength is characterized 
by parameter / such that e£o = f x Uo/2d. For the coupling we consider three regimes: 
(i) a relatively weak field (/ <C 1) where the shape of the quartic potential remains almost 
intact, and the tunneling is still crucially important, (ii) intermediate filed, and (iii) strong 
field / > 1, where the perturbation already considerably changes the low-energy part the 
spectum. 

As the first example, we consider the pulse-driven evolution in a weak magnetic field 
H z = 0.1 T, where the spin dynamics is determined mainly by the SO coupling. The 
dynamics is shown in fig|2J Even a weak field strongly changes the obscrvablcs: in the rich 
transient dynamics it quickly almost equilazes the probabilities to find electron in the right 
and left dot, and the subsequent oscillations of probability around 0.5 value occur. Same 
behavior is demonstrated by the spin components - they weakly oscillate around mean value 
since the electron travels a shorter distance during the density oscillations, however, a similar 
to figJT]spin polarization cr l R (t)/wji(t) is achieved. 

We present driven dynamics in the Hilbcrt space with the coefficients £ n for / = 1/2 in 
figE3 This picture shows only the largest remaining after the end of the pulse components, 
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Fig. 5: Evolution of the probabilities wn(t) (a) and o- z R (t) (b) under the electric field pulse with 
/ = 1/8, / = 1/2, and / = 2 for the initial spin-down state. Circle symbol line corresponds to 
the / = 1/8 regime. The pulse durations T — 25 ps. Similarly to fig[U with the increase in the 
pulse field, the motion becomes less regular, and more high-energy component contribute into the 
dynamics. However, the dynamics is much richer here. For the / = 1/2 case the main underlying 
|£n| 2 are shown in figO 

demonstrating the orbital and spin dynamics, equilibrating the probably density as \^\ 2 
decreases with respect to l^l 2 , and causing the spin evolution by increasing |£i| 2 . 

The obscrvables for the pulses with different 8$ are shown in fig'H] and figO for the 
initial spin-up and spin-down states. Comparison of figO and fig|5] shows that the time- 
dependence of the spin density, and, more important, of the probability is quite different for 
different initial spin states despite the fact that in both cases the dynamics consists of low- 
frequency motion corresponding to the spectrum of the four low-energy states superimposed 
by the set of very-high frequency oscillations due to the involvement of the higher-energy 
states. For the spin-up initial state, the weak and moderate external field just equalizes 
the probability for the electron to occupy the left and the right dots, and then only high- 
frequency oscillations occur. For the initial spin-down state, the dynamics is much richer: the 
low-frequency oscillations with a large amplitude are clearly seen at long times. In both cases 
the amplitude of the fast oscillations increases with the increase in the field. The difference 
between these two regimes is non-trivial and can be understood as follows. The coupled spin- 
charge dynamics is considerably determined by the difference in the energy of the nearest 
spin up and spin down states with opposite parity. The larger is the distance, the less efficient 
is the SO coupling for the coupled dynamics both in the low- and high-frequency domains. 
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For the low-energy spin-up states it is \E& — E±\ = AE g + A z , while for the initial spin- 
down state it is \E a - E 2 \ = AE g - A z . For H z = 2 T, (AE g - A z ) / (AE g + A z ) 0.25, 
making the spin flip at the given spectrum of the pulse, for the latter case easier, and the 
spin dynamics richer. 

To have a better insight into the evolution of the spin density, we illustrate the dynamics 
by a plot of S z (x, t) in ngJBJfor pulse durations T = 12 and T = 25 ps, / = 2 and the initial 
spin-down state. The figure shows that for the shorter T = 12 ps pulse, the S z component 
in the right dot (x = d) is positive at most of the time after the end of the pulse, illustrating 
the controlled spin flip accompanying the electron transfer. By moving along the time axis 
at x = d we follow oscillations of the S z component consistent with the free tunneling on 
the scale of Tr. As for the evolution under the longer pulse T = 25 ps shown in the lower 
panel of figEl the shape of S z (x,t) remains qualitatively the same, however the details of 
the distribution in the right dot are different. Namely, the dominating parts of the profile 
correspond to the negative S z at all times at x « d. This remarkable difference between the 
results for two pulse durations provides a good illustration for the interplay of a driven spin- 
and-chargc motion and the tunneling occurring simultaneously Namely, the stronger and 
the longer the electric pulse is, the more overall changes it provides for the coupled charge 
and spin dynamics. Since the spin evolution is the rotation, a too short (or a too weak) 
electric pulse may " under- rotate" the spin with respect to the desirable state, as we have 
seen earlier for the weak pulses. A too strong (or a too long) electric pulse can "over-rotate" 
the spin, which may pass the desired flipped state in the neighboring dot. This observation 
illustrated in figED suggests to carefully control and adjust both the strength and duration 
of the electric pulses required for the controlled spin manipulation in heterostructures with 
significant SO coupling. 

Conclusions. We have studied the full driven quantum spin and charge dynamics of a 
spin qubit in one-dimensional double single-electron QD with SO coupling. Equations of 
motion in a pulsed electric field have been solved numerically exactly in a finite basis set 
in the regimes of relatively weak, moderate, and relatively strong coupling. The dynamics 
is strongly determined by the pulse amplitude and duration. The results show that the 
time dependence of the spin and charge distributions are strongly irregular with the spin- 
dependent evolution in the charge density pattern. Electron spin flip can be achieved at 
certain pulse durations and amplitudes. These conclusions emphasize the importance of 
the SO coupling which should be taken into account in experimental and technological 
realizations of spin-based qubits. 
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Fig. 6: Evolution of the spin densities under the strong electric field / = 2 pulse for two pulse 
durations T — 12 ps (upper panel) and T = 25 ps (lower panel). The variation in the pulse 
duration leads to the different outcome for the dominating final S z of the transferred electron. 

technical assistance. 
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